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We numerically investigated Mott transitions and mixing-demixing transitions in one-dimensional 
boson-fermion mixtures at a commensurate filling. The mixing-demixing transition occurred in a 
qualitatively similar manner to incommensurate filling cases. We also found the Mott insulator 
phase appeared in both the mixing and the demixing states as the fermion-boson interaction or the 
boson-boson interaction increased. Phase diagrams were obtained in interaction parameter space. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

Ultra-cold alkaline atom gases have opened up new 
possibilities for experimentalists and theorists. The 
trapped atomic gases allow tuning of physical parame- 
ters and provide various kinds of systems from weak to 
strong coupling regimes. For example, commensurate fill- 
ing bosonic systems can be constructed and Mott tran- 
sition has been one of attractive topics studied theoreti- 
cally 0-0 and experimentally 0. 

Easy construction of low-dimensional systems is an- 
other advantage of the trapped atoms. In particular one- 
dimensional boson-fermion mixtures have attracted at- 
tention in terms of mixing-demixing transition, and the- 
oretical studies have been reported 0-01 ■ I* 1 ^ e ^- 
we carried out Monte Carlo simulations to investigate 
the mixing-demixing transition in a system at an incom- 
mensurate filling of fermions and boson, and found that 
the system undergoes the transition from the mixing to 
demixing phases as the boson-fermion interactions in- 
crease. Figure 1 is a schematic phase diagram in the 
Ufb — Ubb plane where Ufb is the fermion-boson interac- 
tion and Ubb the boson-boson interaction. The solid curve 
presents the phase boundary between the mixing and 
demixing phases. An exact solution assures the dashed 
straight line, Ufb — Ubb, is in the mixing phase 0. 

In this paper we examined the Mott transitions and 
the mixing-demixing transitions in systems at a com- 
mensurate filling on a one-dimensional periodic lattice. 
We know single-component commensurate systems be- 
come Mott insulators as inter-particle interactions in- 
crease. In the two-component systems, however, the situ- 
ation is rather complicated partly because the definition 
of the term " commensurability" has some variation. We 
could consider a case where the total filling is commensu- 
rate but the filling of each component is incommensurate 
(Case 1), or a case where the total filling as well as the 
filling of each component are all commensurate (Case 2). 
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FIG. 1: Schematic diagram of mixing and demixing phases of 
incommensurate boson-fermion mixtures. An exact solution 
assures that the dashed line (Ubb = Uf (,)is in the mixing phase. 



Roth and Burnett numerically studied a mixture of 4 
fermions and 4 bosons on an 8-site lattice, i.e. the Case 
2, and showed the presence of mixed Mott insulator and 
demixed Mott insulator at large fermion-boson interac- 
tions and large boson-boson interactions We per- 
formed Monte Carlo simulations to investigate the Mott 
transition and the mixing-demixing transition in the Case 
1 and 2 and see how the phase diagram in Fig. 1 changes 
in the commensurate-filling cases. 



II. MODEL 



We consider a mixture with Nf spinless fermions and 
Nb bosons on an iV-site lattice. To let the system un- 
dergo the Mott transition, we fixed the total number of 
particles N to t — Nf+Nb equal to N. We employed Bose- 
Fermi Hubbard Hamiltonian to describe the interacting 
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particles: 



i a—f,b 



E 



—n b ,i{n bti - 1) + Uf b nf ti n bt ., 



, (1) 



where a a i and CL a i ire respectively creation and annihi- 
lation operators for fermions (a — f) or bosons (a = b) 
on the z-th site, and n a i — Hopping energy and 

chemical potential are denoted respectively by t a and fi a . 
U bb and Uf b are on-site boson-boson and fermion-boson 
repulsive interactions respectively. We set tf = t b = 1 as 
an energy unit. 

We performed Monte Carlo simulations with the world 
line algorithm to measure the system behavior 0, ^| . 
Temperature was fixed to T = 1/15 and the Trotter de- 
composition number L to 100. The number of the sites N 
was set to 14 and the total filling to 1, i.e. N tot = N = 14. 
In the simulation of the Case 1 we chose Nf = 10 and 
N b = 4, and in the simulation of the Case 2 we used 
N f = N b = 7. 

To observe the mixing-demixing transition, we mea- 
sured a correlation function 

C = ((n bl i-i+n bt i +1 )n bl i) - ((«/,»- 1 + + "-/,i+i)«b,i), 

, (2) 

which tends to be smaller, or negatively larger, in mixing 
states than in demixing states. 

To observe the Mott transition, we measured a current- 
current correlation function J to t in the zero-frequency 
limit mill, 



Jtot = lim(jtot(u)jtot(-u)}, 



(3) 



where jtot presents the total current density of fermions 
and bosons. Jtot becomes zero when the system is in the 
Mott insulator phase. 



III. RESULT 
A. Case 1 

At first let us see some limiting cases. In the limit 
of U bb — > oo, i.e. in a hardcore boson case, the mixing- 
demixing transition was not observed and the system al- 
ways stayed in the mixing state. This can be understood 
in the following way. The fermion (boson) gains the ki- 
netic energy by hopping to a boson- (fermion-) occupied 
site, and the gain should be larger in the mixing state 
than in the demixing state, because in the mixing state 
the fermions (bosons) have a lot of chances to hop to 
the site occupied by the boson (fermion). In the limit 
of Uf b — > oo with a finite U bb , on the other hand, the 
system was always in the demixing state. 

What about the Mott transition in these cases? Fig- 
ure 2 shows the behavior of the correlation function J to t 





FIG. 2: The stiffness J to t as a function of Ufb with infinite 
Utb (a) and as a function of Ubb with infinite Ufb (b). 




FIG. 3: (a) The correlation function C as a function of Ufb 
with Ubb = 0. (b) A snapshot of boson configuration (left) 
and fermion configuration (right) at Ufb = 5 and Ubb = 0. X 
is the real space axis and r is the imaginary time axis, rw/l 
shows the number of the bosons (fermions) on each site. 



with infinite boson-boson interactions (a) and that with 
infinite fermion-boson interactions (b). As we can see in 
the figure, the system undergoes the Mott transition at 
U fb - 2 (a), and at U bb - 2.5 (b). 

Next we studied another limit, U bb = 0. The system 
did not show the Mott transition in this case, because 
the bosons can easily overlap each other creating empty 
sites, which makes the system soft. However the sys- 
tem undergoes the mixing-demixing transition. Figure 
3 shows the correlation function C as a function of Uf b . 
The uprising curve of the correlation function indicates 
a transition from the mixing to the demixing states at 
about Uf b = 1.5. The demixing state can be seen in the 
snapshot of Fig. 3(b) where the fermions and the bosons 
stay away from each other. 

In the above limits, the system showed only either of 
the Mott transition or the mixing-demixing transition. 
Now we study the behavior of the system in the param- 
eter region with finite Uf b and finite U bb . In Fig. 4 we 
show the correlation functions C and J, fixing Uf b = 5 
(a) and Uf b = 14 (b). The behavior of C and J indicates 







FIG. 4: The correlation functions C and J tot as a function 
of U bb with (7/6 = 5 (a) and U fb = 14 (b). 



FIG. 5: The correlation functions C and J tot as a function 
of U bb with U fb = 7 (a) and U fb = 14 (b). 



the occurrence of the mixing-demixing transition and of 
the Mott transition both in (a) and (b). The qualita- 
tive difference between (a) and (b) is that the mixing- 
demixing transition occurs almost simultaneously with 
the Mott transition in (a), while the system undergoes 
the Mott transition first and then the mixing-demixing 
transition in (b). 



B. Case 2 

In the Case 2 we have Nt — Nh — 7 on the 14-site 
lattice. As in the Case 1 we measured the correlation 
functions C and J. In the limit of Ubb ^ oo the system 
is in the mixing phase and undergoes the Mott transition 
at around Ufb = 2.5, and in the limit of Ufb — > oo the 
system is in the demixing phase and undergoes the Mott 
transition at around Ubb — 3.5. When Ubb — the system 
changes from the mixing state to the demixing state at 
around Ufb = 1 and never exhibits the Mott transition. 

Figure 5 shows C and J in the intermediate parameter 
region. Their behaviors are quite similar to those in Fig. 
4 of the Case 1. 



C. Phase diagram 

Figure 6 (a) and (b) are the phase diagrams of the Case 
1 and 2 respectively, produced in the Ufb — Ubb plane. It 
is noted that the transition line between the mixing and 
demixing phases is qualitatively equivalent to the one in 
the case of incommensurate filling as in Fig. 1. The phase 
diagrams are also consistent with the result in Ref. Q 
which proved that the Ufb = Ubb line stays in the mixing 
phase. 

We do not find any essential difference between the two 



diagrams (a) and (b), indicating that the phase diagram 
is determined by the commensurability of the total fill- 
ing and not very sensitive to whether the filling of each 
component is commensurate or not. 

Figure 7 shows the snapshots of typical configurations 
of the Case-1 particles in the mixed Mott insulator (a) 
and in the demixed Mott insulator (b). The particles fill 
every site uniformly and make the system insulating both 
in (a) and (b), but the fcrmions and bosons are mixed in 
(a) while separated in (b). 



IV. CONCLUSION 

We investigated one-dimensional boson-fermion mix- 
tures with one particle per site. Focuses were placed 
on two cases. In the Case 1 the fermion filling and bo- 
son filling are both incommensurate but the total filling 
is commensurate, and in the Case 2 the fermion filling 
and boson filling are both commensurate. The Monte 
Carlo simulations revealed the occurrence of the mixing- 
demixing transition and the Mott transition. The system 
showed the mixing-demixing transition at a certain value 
of Ufb and the transition point shifted to the larger Ufb 
side as Ubb increases, which was qualitatively similar to 
the incommensurate filling case [ij. We found the Mott 
transition occurred in both mixing and demixing phases 
at a sufficiently large interaction. The phase diagrams of 
the Case 1 and 2 were obtained and clearly exhibited the 
presence of the mixed Mott insulator and the demixed 
Mott insulator. 
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FIG. 6: Phase diagrams of the Case 1 (a) and Case 2 (b). 
Ufb is infinite on the right axis and Ubb is infinite on the upper 
axis. The filled diamonds present the Mott transition points 
and the empty squares the mixing-demixing transition points. 
The solid lines are just the guides for eyes. MI denotes the 
Mott insulator. 
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FIG. 7: Snapshots of the particles in the Case 1. (a): The 
mixed Mott insulator with Ufb = 7 and Ubb = 6. (b): The 
demixed Mott insulator with Ufb = 14 and Ubb = 4. The 
upper figure shows the boson configuration and the lower the 
fermion configuration. 



